Wigner Crystallization in Rapidly Rotating 2D Dipolar Fermi Gases 
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We study the competition between the Wigner crystal and the Laughlin liquid states in an ul- 
tracold quasi two-dimensional rapidly rotating polarized fermionic dipolar gas, and find that the 
Wigner crystal has a lower energy below a critical filling factor. We examine the quantum crystal to 
liquid transition for different confinements in the third direction. Our analysis of the phonon spectra 
of the Wigner crystal with the account of phonon-phonon interactions also shows the stability of 
the Wigner crystal for sufficiently low filling factors (y < 1/7). 



Recently a remarkable progress has been made in ex- 
perimental studies of strongly correlated systems of ul- 
tra cold gases. A seminal breakthrough was the obser- 
vation of Mott insulator-superfluid transitions in Bose 
gases in optical lattices [l| , followed by studies of low di- 
mensional gases, fermionic superfluidity or ultracold dis- 
ordered gases (for recent reviews of this very rapidly de- 
veloping area see 0,9). A particularly fascinating route 
toward creation of strongly correlated states is the one 
that uses rapidly rotating gases. The rotation is formally 
equivalent to a magnetic field (in the rotational frame), 
which reorganizes free-particle states into discrete highly 
degenerate Landau levels. As a result, the properties 
of the system become very sensitive to interparticle in- 
teractions, in complete analogy to the Fractional Quan- 
tum Hall Effect (FQHE) for electrons 0], that exhibits 
a variety of strongly correlated states, among which the 
Laughlin liquid [5(, and the Wigner crystal |6[ are prob- 
ably the most famous. This analogy has been pointed 
out in the context of ultracold atomic gases with short 
range interactions in Refs. 0|, but an experimental ob- 
servation of the FQHE in this case is difficult, partially 
due to smallncss of the energy gap separating a Laughlin 
state from the so called quasi-hole excitations. It was 
recently argued [f| that these problems may be overcome 
in rotating quasi-2D polarized dipolar gases, in which 
the interparticle interaction has a long-range part that 
puts dipolar systems somewhat in between atomic sys- 
tems with short-range interactions and electron systems 
with Coulomb interactions. As a result, the Laughlin 
state for the filling factor v = 1/3 [32| was predicted to 
be incompressible with gapped excitations [8[ . 

Several other unexpected predictions have been al- 
ready made for rotating dipolar bosonic gases. The ro- 
tating dipolar condensates in the mean-field regime ex- 
hibit novel forms of vortex lattices: square, "stripe crys- 



tal", and "bubble crystal" lattices 0, [l(|. On the other 
hand, Rezayi et al. (TH have recently shown that a 
small amount of dipole-dipole interactions stabilizes the 
so called bosonic Rezayi-Read strongly correlated state 
at v = 3/2 whose excitations are both fractional, and 
non-Abelian. 

The dipolar interactions could also lead to a crystal 
ground state in a quasi 2D system similar to the electron 
Wigner crystal. In the case of an electron gas, both in the 
presence and in the absence of a magnetic field, Wigner 
crystal becomes energetically favorable and stable for low 
filling factors: this is due to the fact that at low densi- 
ties Coulomb interactions (~ 1/R, where R is a typical 
leng th scale) win over the kinetic energy (~ 1/-R 2 ) (cf. 
[l2j). For non-rotating dipolar gases (where interactions 
scale as l/R 3 ), it was argued [8( that the ground state 
has a crystal order at high densities, as shown in Refs. 
[l3j ]. Remarkably, in the case of a rotating dipolar gas, 
the Wigner crystal phase is expected at low densities [f| . 
Therefore, the dipolar gases, being to some extend similar 
to electrons, present nevertheless a novel class of physical 
systems exhibiting quantum liquid to crystal transition. 
This opens an unprecedented possibility for the observa- 
tion and detailed studies of this transition in clean and 
well-controlled conditions. It is especially appealing in 
view of the recent experimental realization of the dipolar 
Bose condensate of Chromium [3, EH , a nd p rogress in 
trapping and cooling of dipolar molecules [1 61 ] - 

In this paper we investigate the existence of a Wigner 
crystal phase in a rapidly rotating gas of polarized dipo- 
lar fermions by comparing it with the competing liquid 
states. We then examine the stability of the Wigner crys- 
tal by incorporating phonon-phonon interactions, and 
identify the border of stability with the appearance of 
purely imaginary phonon frequencies. It appears that 
the Wigner crystal is the stable ground state of the sys- 
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terns for low filling factors, typically for v < 1/7. 

We consider a gas of fermionic polarized dipolar par- 
ticles in a rotating cylindrical trap, tightly confined and 
polarized along the axis of rotation. For a sufficiently 
strong confinement in the axial direction, all particles 
are in the ground state of the axial motion and we can 
write the many-body wave function in the form 

iMfa.f,}) = lM{r*}) {lV^)- N, \M-Y.^il 21 ^ 

where = (x, y) and & are radial and axial coordinates 
of the particles, respectively, I is the extension in the 
axial direction, N is the total number of particles, and 
ip2D is the wave function of the system in the xy-planc. 
The confinement in the axial direction also modifies the 
interparticle interaction in the 2Z)-plane. For two dipolcs 
separated by a distance r in the 2D-plane, the effective 
interaction potential is 
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For r 3> I the potential is i>2D ~ d 2 /r 3 , while for r < I 
the effective potential increases logarithmically, V2U ~ 
(d 2 /l 3 ) y/2/nhi(l/r). The wave function ijj2D obeys the 
Schrodinger equation with the effective 2D Hamiltonian 
that in the rotational frame reads 
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where Q is the rotational frequency, u>± is the radial 
trap frequency, Vd = J2i<j v ?d (|r, — Tj\) describes the 
dipolc-dipole interparticle interaction, and L z is the z- 
componcnt of the total angular momentum. The Hamil- 
tonian ([1]) can be rewritten in the form 

H = V^-HRV,- A,) 2 - (Q-u±)L z + Vd (2) 
^— ' 2m 



with Aj = mu±_e z x and, therefore, it formally de- 
scribes a system of charged particles in a constant mag- 
netic field with the cyclotron frequency uj c = 2ujj_. For 
non-interacting particles (Vd = 0) in the regime of the 
critical rotation (uj± = Q), the properties of the Hamilto- 
nian |(5J) arc well known; it's spectrum consists of highly 
degenerate levels E n = hui c (n + 1/2) called Landau lev- 
els. In the following, we restrict oursclf to a system of 
particles occupying only the lowest Landau level. 

The wave function of a non-correlated Wigner crystal 
on the lowest Landau level is [T3| : 



*c({z,}) ~-4jJexp 



where Zj = x% + iyi is the complex representation of the 
2D vector r,, Ri the (complex) lattice site, Iq = \J h/muj c 



the magnetic length, and A the antisymmctrization op- 
erator, which can be omitted for sufficiently low fill- 
ing factors v = 2-nl\n. One can check that similar to 
the case of classical dipoles, the energy is minimal for 
a triangular lattice with particles centered at positions 
Ri = Zibi + l%\>2-, where l\ % 2 are integers, bi = a (0,1), 
b-2 = a (s/3, l) /2, and a is the lattice constant deter- 
mined by the density n of the gas, a 2 = 2/v / 3fi = 
47tZq/v / 3^- The energy of the crystal state equals 

U c = (*c|Vb|*c) ~ d 2 n 3/2 (0.2823 + 0.2146/3 
+ 0.3388/3 2 + 0.7456/3 3 + 2.06 76/3 4 + ...), 

where (3 = 7rn,(2Z 2 — I 2 ) and the first term corresponds 
to the energy of classical point-like dipoles, E c \ = 
Z t<j rf 2 /l^ibi + hh 2 \ 3 = 5.513d 2 /a 3 . 

The energy of the Laughlin state (for a small number 
of particles interacting via dipolar forces, recent compu- 
tations show a remarkable ove rlap of the exact ground 
state with the Laughlin state [18|) characterized by an 
odd integer M = 1/u is 



rg{r)v2D {r)dr, 



where g(r) is the pair correlation function. We calculated 
this functions and the energies for odd M from 1 to 19 
for a gas of 512 particle using a Monte Carlo method. 

We now compare the energies of the Wigner crystal 
and of the Laughlin liquid for different filling factors v 
and extensions I in the axial direction. The results are 
presented in Fig. [T] They show that below some critical 
value u c , which depends on /, the Wigner crystal has a 
lower energy than the Laughlin liquid and, therefore, for 
v < v c the ground state is expected to be a crystal. Note 
that in obtaining the critical filling factor we compare the 
energies of two simplest trail wave functions for liquid 
and crystal states. A better estimate could follow from 
considering the wave functions of quantum Hall liquids 
of composite fermions [l^] with the filling factors closer 
to the critical one than primary Laughins 1/M, and of 
the correlated Wigner crystals [20j. We, however, do not 
expect a significant change of our result because for a 
dilute system the two considered wave functions already 
catch the most important features of the interparticle 
interaction. 

Let us now approach the liquid-crystal transition start- 
ing from the crystal phase. For this purpose we consider 
the phonon spectrum with taking into account phonon- 
phonon interactions (anharmonicity effects). The ap- 
pearance of purely imaginary phonon frequencies will in- 
dicate the instability of the crystal and, therefore, the 
transition to a liquid state. 

In the harmonic approximation the phonon eigenfre- 
quencies can be obtained from the dynamic equations 
for displacements u a \ of particles along the a-axis from 
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Figure 1: Energy per particle for the Wigner crystal (dotted 
line) and for the Laughlin state (solid line) as a function of 
the filling factor for I — 0. The insert shows the critical filling 
factor as a function of the extension in the z-direction. 



their equilibrium positions Ri in the lattice 21 1 
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where sv = d 2 Uc / dR a \dRj3\i is the dynamical ma- 
trix, a, (3 = x,y and e a p is the antisymmetric tensor, 
e xy = 1. The last term in Eq. ([3]) corresponds to the 
Coriolis (Lorentz) force due to the rotation (magnetic 
field). 

Without rotation, lj c = 0, Eq. [3] in the quasimomcn- 
tum representation reads: 

uj 2 u a (k) = ^2 F a p(k)up(k), 



where mF aa i (k) 
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Fourier transform of $ 
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and the eigenvalues wf(k), 
s = T, L of F aa i (k) determines the frequencies of the 
transversal (T) and longitudinal (L) phonons. The 
phonon frequencies for I = are shown in Fig. [2] They 

lot ~ (3 / \fS)iOdak and 



are linear in k for k <C a -1 

ujl ~ VTLujt, where u>d = ^d 2 /ma 5 sets the typical 
value for the frequency of phonons. 

For a rota ting crystal, lu c ^ 0, the phonon spectrum 
becomes [2lL l22l 
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In the case lo c 3> oj s , one has oj+ m 
ljlut/uc- 

Higher orders (anharmonic) terms in the expansion of 
the energy of the crystal with respect to the displace- 
ments of particles from their equilibrium positions re- 
sult in phonon-phonon interactions and, therefore, in 
the renormalization of the phonon frequencies. At 
a given quasimomentum k, the renormalizcd frequen- 
cies correspond to the poles of the Fourier transform 
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Figure 2: Energy of the transverse ujt (lower surface) and 
longitudinal ujl phonons. 

Figure 3: One-loop diagrams for the phonon self-energy. 



G a p(uj, k) of the phonon Green function G a p(t, 1 — 1') = 
i [9{t){ Ucd {t)up V {Q)) + 9(-t)(u 0V (O)u a i{t))] with 9{t) be- 
ing the step function (for more details on the usage 
of Green functions for phonons in crystals see e.g. 
[23I]). The Green function in the harmonic approxi- 



mation is G^l 1 (uj,k) 



ie a f3Ld c oj}m/h, where = ea ! e y ^' is the projector to 

the eigenmode s with the polarization e^ s \ and the poles 
are at lo = uj±(k). In the presence of phonon-phonon in- 
teractions, the Green function obeys the Dyson equation 



:^ 2 (k) 



G- 1 (o J ,k) = G^- 1 (o J ,k)-E a0 (w,k), 

where the phonon self-energy function E incorporates all 
the effects of phonon-phonon interactions. 

We calculate the phonon self-energy in the one-loop 
approximation that has proven its reliability in the case 
of an electron Wigner crystal in a strong magnetic field 



24, 25, 26]. This corresponds to taking into account only 
the fourth $( 4 ) - d 4 U G /dR i and the square of the third 
<j>( 3 ) r^, d 3 Uc/dR 3 order anharmonic terms, see Fig [3] 
We then solve the Dyson equation numerically by suc- 
cessive iterations until we either obtain a self-consistent 
solution or the iteration breaks down due to the appear- 
ance of purely imaginary phonon frequencies. The lat- 
ter indicates the instability of the lattice. The results 
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of our calculations can be summarized as follows. At 
very low filling factors the phonon-phonon interactions 
do not play any significant role and the system behaves 
harmonically. However, the effects of anharmonicity be- 
come progressively important with increasing the filling 
factor, and, finally, at some critical filling factor v c we 
find purely imaginary frequencies and, hence, the crystal 
ceases to exist. We always observe this phonon instabil- 
ity for k — > 0. This indicates the breakdown of the crys- 
talline order and rules out any possible structural phase 
transition. Therefore, the value v c sets the upper bound 
for the stability of the crystal. The critical filling factor 
v c depends on the confinement in axial direction as well 
as on the ratio of u>d and lu c that measures the strength of 
the dipole-dipole interaction relative to the Landau level 
spacing. This dependence for / = can be approximated 
as v^ 1 = 4.33exp(-0.0021w c /cj d ) + 5.77. Note that for 
bj c 3> LOd the critical filling factor v c becomes insensitive 
to the strength of the interparticle interaction. In this 
limit the Wigner crystal is stable for v < v c = 0.174, 
which is in a good agreement with our previous energetic 
consideration for the tight confinement along the z-axis, 
I < l . 

We also calculate the Lindemann parameter 7 - the 
ratio of the average displacement of a particle in the lat- 
tice from its equilibrium position to the lattice spacing 
a, 7 = yj (u 2 )/a. For lo c 3> o->d we find 7 = 0.28 that is 
within the range of values of the Lindemann parameter 
for various 2D crystals. 

As we see, the crucial parameter that controls the 
ground state of the system and, therefore, drives the 
quantum melting transition is the filling factor. The way 
of manipulating the filling factor in experiments with ro- 
tating gases depends on an experimental setup. In the 
case of a critical rotation with an extra (quartic) confine- 
ment [llj this can be achieved by changing the number of 
particles (the size of the system is fixed by an extra con- 
finement). For a purely harmonic confinement and "un- 
der critical" rotations [2^, [29| , the filling factor can be 
changed by varying the difference f2 — In this setup, 
the size of the system and, therefore, the filling factor is 
determined by a competition between the interparticle in- 
teraction Vb and the "tilting" term (fi — ui±) L z (see Eq. 
[2]) . In both cases, the appearance of a crystal order could 
be detected by studying the shot noise correlations using 
the Hanbury Brown and Twiss effect (30j in a similar way 
as it was used to observe the Mott insulator-superfluid 
transition in a lattice Bose gas [3lj | . 

In conclusion, we show the existence of the Wigner 
crystal state for low filling factors in a rapid rotating 2D 
dipolar Fermi gas. The critical filling factor, at which the 
quantum melting of the crystal into a quantum Hall liq- 
uid takes place, depends on the confinement in the axial 
direction and on the strength of the dipole-dipole inter- 
action, but typically one expects a liquid ground state 
for v > 1/7. 
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